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Based on spatial-temporal symmetry breaking mechanism, we propose a novel scheme for terahertz 
(THz) wave generation from hyper-Raman lines associated with the 0th harmonic (a particular even 
harmonic) in a two-level quantum system driven by two-color laser helds. With the help of analysis of 
quasi-energy, the frequency of THz wave can be tuned by changing the field amplitude of the driving 
laser. By optimizing the parameters of the laser fields, we are able to obtain arbitrary frequency 
radiation in the THz regime with appreciable strength (as strong as the typical harmonics). Our 
proposal can be realized in experiment in view of the recent experimental progress of even-harmonics 
generation by two-color laser fields. 

PACS numbers: 42.65.Ky, 42.50.Hz, 78.20. Bh 



Introduction Terahertz (THz) radiation, electromag- 
netic radiation with typical frequency from 0.1 THz to 10 
THz-lying in the spectrum gap between the infrared and 
microwaves, has varieties of applications in information 
and communication technology, biology and medical sci- 
ences, homeland security, and global environmental mon- 
itoring etcji. The generation of THz wave is a challenge 
problem due to the lack of appropriate materials with 
small bandgaps in the usual optical approach. Great ef- 
fort has been made to the design of THz sourcesSrJ^. 
THz wave can be obtained by both electronic and opti- 
cal methods. The uni-travelling-carrier photodiode^ and 
the quantum cascade laser—"— are two examples. The 
electronic approach is limited in the low frequency end 
of the THz regime, and the optical approach usually fo- 
cuses on the cases of small energy gap. Technological in- 
novation in photonics and nanotechnology has provided 
us more new ways for THz radiation generation. Re- 
cently an interesting approach based on semiconductor 
nanostructure driven by acoustic wave was suggestedi^, 
and electrically pumped photonic-crystal THz laser was 
developedii. THz wave generation by up conversion 
method through high-order harmonics generation(HHG) 
in semiconductor quantum dot(QD) was also suggested 
in our previous worki^. 

One challenge of THz wave generation in optical ap- 
proach in usual system is to generate low frequency ra- 
diation from a quantum system with a large energy gap. 
Usually, in the emission spectrum of a quantum system, 
there are harmonics, as well as the associated hyper- 
Raman lines. Hyper-Raman lines are caused by the tran- 
sitions between the dressed bound states. Nth hyper- 
Raman lines are defined as those associated with the nth 
harmonics (they appear with the corresponding harmon- 
ics simultaneously due to the same symmetry reason, and 
locate near uloq, where loq is the fundamental frequency 
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of the incident laser Here is one observation that the 
shift of the hyper- Raman line to the associated harmonic 
could be small and could be tuned by the external field, 
even though the fundamental frequency is large. Then 
one may use the 0th hyper-Raman line (that associated 
with the 0th harmonic) to generate the low-frequency 
(THz) wave. There have been some studies on the har- 
monics and hyper-Raman lines. In many cases, there 
are only odd harmonics due to the particular symmetry 
and the associated hyper- Raman line is of high frequency 
and/or weak intensity. Therefore one has to solve these 
two problems to generate THz wave by using the hyper- 
Raman lines in the nonlinear optical processes. 

In this article, we adopt an effective optical method to 
obtain THz wave in typical two-level quantum systems, 
which solves the two problems mentioned above. Based 
on the generalized spatial-temporal symmetry principle 
we developed recently^S, we propose a novel mechanism 
to obtain the 0th hyper-Raman lines in THz regime by 
introducing additional laser with frequency 2fcaJo (k = 
1,2,3, ...). Using analysis of quasienergy and optimiza- 
tion of laser fields, we are able to obtain considerably 
intense THz radiation with desired frequency. With the 
help of our optimization method, it is quite likely that 
our proposal can be realized in experiment considering 
the recent experimental progress of even-harmonics gen- 
eration by two-color laser. 

Theoretical formulism Our two-level quantum sys- 
tem driven by laser fields [see the schematic diagram Fig. 
1(a)] is described by the Hamiltonian 

2 

i/ = ^i?,|*)(z| + G(0(|l)(2| + |2)(l|), (1) 

i 

where G{t) ~ F{t)e ■ /xi2, is the Rabi frequency caused 
by periodic laser field E = F{t)e, F{t) = F{t + T), e a 
unit vector. /ii2 = (l|er|2) is the dipole between state 
|1) and state |2). The energy spacing between the two 
states [with energies Ei{i = 1, 2)] is set as AE = Ei—E2- 
The dynamics of our system is described by the equa- 
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tion of motion of the density matri»21 
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where H is the Hamiltonian, the last term describes 
possible dissipative effects (such as spontaneous phonon 
emission) and we set ?i = 1 in the following. We nu- 
merically calculate the photon emission spectra of the 
two-level quantum systems by solving the density ma- 
trix in Eq.(I2|) through Runge-Kutta method with the 
time step of 0.002/a;o, total steps of 1200000 and the 
electron initially setting at the lower level. The average 
dipole can be calculated as D(i) = J2ij f^ijPij{t)- We use 
Fourier transformation to obtain the emission spectrum 
S{i') = I Jdtex.p{—ii't)'D{t)\^. In the calculation we set 
hujo (cjo the driving field frequency) as the unit of energy. 
We choose ujq = lOOTHz as an example, yet in general 
one can also use other frequency (much larger than THz) 
laser as will be discussed later. The two-level quantum 
systems can be realized in many systems, such as atoms, 
molecules, and semiconductor quantum dots. Here we 
use the typical parameters of a quantum dot: the energy 
spacing AE = IOHluq, the dipole moment /ii2 = 0.5 e-nm, 
phonon emission coefficient Tij = 3.4 x 10~'^a;o. 

We explore the emission spectra using above numer- 
ical approach and the theory on the selection rule in 
high-order harmonic generation based on the generalized 
spatial-temporal symmetry. As described in our previous 
work2£, for a quantum system described by Hamiltonian 
(1), if there exists one symmetric operation Q, which is 
the time shift 9 : t ^ t + T/2 (T = 27r/wo) combined 
with another operation in spatial/spectrum domain, 
i.e. Q = ■ 9 , such that the initial condition and the 
Hamiltonian (up to a sign) arc invariant, and the dipole 
operator P has a definite parity, then the emission spec- 
trum contains no odd/even component if operator P is 
even/odd. For our two-level quantum system driven by 
a monochromatic laser, i.e. E{t) = Ficos{uJot), there ex- 
ists one symmetric operation Qi, which is the time shift 
9 combined with spatial operation ci ^ — ci (here 
and in the following Cj, j = 1,2, refers to the annihilation 
operator for state |j)), and the Hamiltonian is invariant 
and the dipole operator P is odd. Then odd harmonics 
alone are generated because of the spatial-temporal sym- 
metry. Other types of spatial-temporal symmetry may 
lead to many interesting emission patterns^S. 

Here we give some analysis on the hyper-Raman lines 
and their relation to the harmonic components. In our 
systems driven by a periodic field, the quasiencrgy state 
has the form \i/ja{t)) = e~"^°*|(/)a(t)), where the Floquct 
state \(j)a{t)) ~ \(f>a{t + T)) can be written in the form 
\(j)a{t)) = J2m e"™"''*!^)- We consider a state |?/'(t)) = 
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FIG. 1: (a) The schematic diagram of our two-level sys- 
tem driven by an incident laser. The emission spectrum 
contains the higher-order harmonics and the accompanied 
hyper-Raman lines; (b)-(d) Emission spectra in the presence 
of various driving fields F{t). (b) F{t) = Fi cos{u}ot); (c) 
F{t) = F2cos{2ujot); (d) F{t) = Ficos{ujot) + F2 cos(2LJot); 
(e) F{t) = Fi cos(cjoi)+-f2 cos(3ajoi); (f) F{t) = Fi cos{ujot) + 
F2Cos{iLoot). Fi = 4.4 X 10^V/m,F2 = 2.2 x lOV/m. 



If the system posses a symmetry generated by the oper- 
ator Q = Vl-9 and the Floquet state |(/'a(i)) has a definite 
parity under Q, i.e., Q\4>a{t)) — ±|(/)Q(t)}, then we have 
= ±(-l)™|0™). Then we have (Cl-P|<^^) = 
(a, /3 = 1, 2) for n~m even/odd number and states a,l3 of 
same/different parity. Quite often the harmonics and the 
associated hyper-Raman lines appear simultaneously due 
to the same symmetry properties of the Floquet states. 
In the case with a monochromatic laser, we have odd har- 
monics and the associated hyper-Raman lines. While in 
other cases with symmetry broken, more harmonics and 
the associated hyper-Raman lines are generated. 

Emission patterns Let's first look at the emission 
spectrum of a system driven by a monochromatic inci- 
dent laser with frequency ujq. As seen in Fig. 1(b), there 
arc odd harmonics as well as the associated hyper- Raman 
lines due to the spatial-temporal symmetry properties of 
the Floquet states as analyzed above. It is natural that, 
for the incident laser with frequency 2a;o, the emission 
spectrum contains components with frequencies of 2wo, 
QuJo and 10a;o...(odd orders of incident frequency 20^0) 
as shown in Fig. 1(c). One notices that there is no 0th 
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hyper-Raman line in this case. If we use two-color lasers 
with frequencies loq and 2ujq, interesting phenomena ap- 
pear. As seen from Fig. 1(d), even harmonics and the 
associated hyper- Raman lines are generated by introduc- 
ing the second laser field. Especially a low frequency 
radiation, 0th hyper-Raman radiation is generated. We 
would like to point out that the emission spectrum for the 
case with driving field F{t) = Fi cos(a;ot) + F2 cos(2ajoi) 
is not the supposition of those driven by Fi coa{ujot) and 
F2 cos{2uJot). For example, the harmonic 4wo in Fig. 1(d) 
neither appears in Fig. 1(b) nor in Fig. 1(c). It is also not 
the frequency summation or difference of the harmonics 
of systems driven by monochromatic incident laser with 
frequency ujq and 2wo- In fact, the appearance of all even 
components (and the associated hyper-Raman lines) in 
Fig. 1(d) is the consequence of symmetry breaking as 
discussed above. 

According to our theory, the symmetry of the quan- 
tum system generated by Qi is broken by introducing 
the second laser with frequency 2fcwo {k — 1,2, 3...), and 
it is not broken by introducing the second laser with fre- 
quency {2k — l)wo {k = 1,2,3...). These predictions are 
verified by our numerical results shown in Fig. 1(e) (the 
second laser of frequency 3a;o) and Fig. 1(f) (the second 
laser of frequency 4a;o). It is clear that the 0th hyper- 
Raman line does not appear in the cases with symmetry 
(see Figs. 1(b) (c) (e)). Interestingly, the second laser field 
with frequency Aluq leads to the appearance of 2nth har- 
monics (even for n ~ I) and the associated hyper-Raman 
line (in particular the 0th hyper- Raman line) as predicted 
by our theory, since the second laser with frequency 4a;o 
breaks the spatial-temporal symmetry generated by Qi. 
Using this mechanism, we can explain the experimental 
results of observing even harmonics in helium or plasma 
plumes (containing nanoparticles, carbon nanotubes, etc) 
driven by a two-color laser,— where the second-harmonic 
driving term, despite being very small, breaks the sym- 
metry and thus allows additional strong even harmonic 
components. Our theory also naturally explains the gen- 
eration of even harmonics and associated hyper-Raman 
lines by breaking the spatial symmetry^. 

As seen from equation (3), the frequency of hyper- 
Raman line is determined by the quasienergy, which is 
tunable. One may ask whether one can use the hyper- 
Raman line associated with 1th harmonic to generate 
the low frequcncy/THz wave? Actually, as seen from 
Fig. 2(a) the intensity of the hyper-Raman line asso- 
ciated with the 1th harmonic decreases dramatically as 
the frequency decreasing. To make this point clearer, we 
compare the low frequency components associated the 
1th harmonic and 0th harmonic for the case with two- 
color laser shown in Fig. 2(b). One sees that in the ev- 
ery low frequency regime, the intensity of the 0th hyper- 
Raman line (associated with the 0th harmonic) is much 
larger than that of the 1th hyper-Raman line. Therefore 
one should use the 0th hyper-Raman line to generate the 
THz radiation effectively. 

Tuning of the frequency and intensity of THz 
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FIG. 2: (a) The intensity of hyper-Raman line associated 
with 1th harmonic for F{t) — Fi cos{u;ot), F\ is related to v, 
the frequency of the hyper-Raman line; (b) The intensities 
of hyper-Raman lines associated with 0th (peak 1) and 1th 
(peak 2) harmonics for F{t) = F\cos{u}ot) + F2Cos{2u)Qt), 
Fi = 4.4 X l(fV/m, F2 is related to u. 




FIG. 3: The quasienergy difference (sohd line) between two 
quasi-eigenstates and the frequency (dots) of the 0th hyper- 
Raman line versus the magnitude of external field -Fi. F2 = 
2.2 X lOV/m. 



wave Now we discuss the optimization scheme of THz 
wave generation with driving field F{t) = Fi cos(a;oO + 
-^2 cos(2u;oi) in our proposal. From equation (3), we see 
that the difference of the quasienergy is related to the 
shift of the hyper-Raman line form the corresponding 
harmonics, or the frequency of the 0th hyper- Raman line. 
This relation is verified by our numerical calculation as 
shown in Fig. 3. Moreover from Fig. 3, one sees that 
the 0th hyper-Raman line with very small frequency can 
be obtained in some parameter regimes, which can be 
used for THz wave generation. By applying this rela- 
tion A = I El — £2! = TT'V [v the frequency of the 0th 
hyper- Raman line), any designated frequency THz wave 
can be obtained by tuning the magnitudes of the two ex- 
ternal fields Fi and F2- For example, THz radiation of 
2THz (0.02ajo) can be obtained by tuning the external 
field as shown in Fig. 3 [intersection points of the solid 
straight-line {v = 0.02a;o) and the line for quasienergy]. 
Fig. 4(a) shows the values for Fi and F2 under which 
the low frequency is = 0.02a;o. There are disconnected 
parameter regimes of driving fields for generation of THz 
radiation. Furthermore, the intensity of the THz wave 
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FIG. 4: (a) The parameters of Fi and F2 for the 0th hyper- 
Raman hne with frequency 2THz= 0.02a;o- The triangle indi- 
cates the optimal parameter for the maximal intensity of 0th 
hyper-Raman line, (b) The corresponding emission spectrum 
for the case with optimized external field. 



can be optimized. One can find the proper parameters 
for the maximal intensity with desired frequency [the tri- 
angle in Fig. 4(a)]. The corresponding emission spec- 
trum is shown in Fig. 4(b). It is seen that the intensity 
of THz radiation has been increased 200 times after op- 
timization and we are able to obtain appreciably intense 
THz wave (as strong as typical harmonics). Here we es- 
timate the emission power based on a system of arrays 
of semiconductor quantum quantum dots. The emission 
power of each dot is P ~ l/xispi^^'/Sc^ - lO'^^W. There 
are about iV ~ 10^ dots in the regime of size of wave- 
length which emit wave in phase. Therefore, the total 
emission power from a sample of submillimeter size is 
around N^P ^ fiWM 

In our approach, the typical driving field intensity 
is in the order of 10^^ ^ lO^^W/cm^ , which is lower 
than the typical driving field intensity (in the order of 
lO^^W/cm'^ lO^^W/cm'^) used in other methods for 
THz generation and/or HHG by two-color laser— 1^—. 
One should also use a long driving laser pulse [much 
longer than its period {2tt/ujo)] as that used in Ref. [2^ . 
The emitted hyper-Raman line frequency can not only 
be tuned by the driving field intensity as seen in Fig. 
3, but also by the frequency of the incident laser. If 
we use the incident laser with frequency in the regime 
ujQ ^ hOTHz — 500THz, we may obtain the 0th hyper- 
Raman line with frequency -fqOJo ^ ITHz — lOTHz. 



In our method, there is another tunable parameter, the 
phase difference between two incident lasers. The main 
physical picture remains the same for different phases. 

Unlike previous studies on non-linear driving two-level 
systems^^, the hyper-Raman lines and low frequency 
generatio n^^i^^" — , our scheme of generating low fre- 
quency (THz) wave is based on the symmetry principle in 
typical two-level quantum systems of large energy spac- 
ing (in the order of eV). It can be easily obtained from 
natural systems (such as atoms, molecules, and semicon- 
ductors) and artificial structures (such as quantum dots) . 
The large energy spacing makes it robust against ther- 
mal fluctuation and also insensitive to the initial stateSS. 
Our mechanism based on the symmetry principle is dif- 
ferent from those (based on four-wave mixing^ or tran- 
sit current^) used before. Our non-perturbative anal- 
ysis and calculation show that one only needs to tune 
the intensity of the two-color laser for typical two-level 
quantum systems, which is very convenient in experi- 
ments. We may be able to obtain tunable THz radiation 
due to the availability of stable, tunable driving sources 
above the THz regimeii^. Moreover, people have already 
observed even harmonics in experiments with two-color 
laser— Therefore it is quite likely that the 0th hyper- 
Raman line (in THz regime) associated with 0th har- 
monic can be obtained in experiment if our optimization 
method is used. 

Summary A novel mechanism based on spatial- 
temporal symmetry breaking is proposed to obtain THz 
wave from Otli hyper-Raman line in two-color pumped 
two-level quantum system. Quasienergy is calculated to 
determine the parameters of the incident laser to obtain 
radiation of arbitrary frequency from O.lTHz to lOTHz. 
Upon optimization of the driving fields, we are able to 
obtain THz wave with desired frequency and appreciable 
intensity (as large as that of the typical harmonics). 
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